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I. INTRODUCTION 



The study of the implications of quantum mechanics applied to multipartite correlations has brought many apparent 
paradoxes [1, 2] and later many useful applications in information processing [3], communication [4], cryptography [5] 
and metrology [6]. One of the simplest systems that is able to show a lot of the complexities and beautiful subtleties 
of quantum mechanics is a pair of two level systems (or qubits) . As simple as it may look, many important questions 
regarding this simple bipartite system have not been or cannot be answered [7] . While the situation for the description 
of the correlations of pure quantum states is reasonably well understood, extending the treatment to encompass mixed 
states, in which quantum and classical correlations both play roles, has led to a veritable menagerie of metrics intended 
to gauge the "quantumness" of a state. In order to clarify things, in this paper we aim to provide a compendium of 
the most relevant properties of certain type of two qubits systems known as X states. These states are ubiquitous in 
the literature as they generalize many important classes of mixed quantum states such as maximally entangled states 
(like the singlet state or the Bell states), partially entangled and quantum correlated states (like the Werner states), 
the maximally entangled mixed states [8] as well as non entangled, non quantum-correlated states. The states that 
concern us here receive their name for the form of their density matrix: 
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Px = z* c ' (1) 
V w* d J 

in the basis 1 00) , 1 01} , 1 10) , 1 1 1) (\a/3) = \a) A (g> \/3) B ). They have been generated in different physical systems. Two 
of such systems are the polarizations of a pair of photons generated by a non-linear crystal [9] and the electronic 
levels of a pair of cold ions in a trap [10]. In both cases pure states are generated. Nevertheless, one can imagine that 
by a desired process[ll] (such as having the photons pass through a decoherer) or an undesired one [12] (like stray 
magnetic fields that randomly shift the local energy levels of the cold ions) only coherences between different basis 
elements are preserved. 

For instance it might happen that a certain pure two qubit state vector \ip) — a|00) + /?|01) + 7|10) + S\U) is 
prepared. These qubits can be encoded in the levels of a pair of trapped cold ions. In this case, the levels will be 
subject to the Zeeman effect induced by stray magnetic fields and they will suffer a shift 

AE = n B gjmjB(t) (2) 

where /is is the Bohr magneton, gj is the gyromagnetic factor, mj is the magnetic quantum number of the level and 
B(t) is the fluctuating magnetic field. Under these circumstances the relative phase between the excited (|1)) and 
ground (|0)) will be given by: 



Hi) = f { 9 fmf- g fm^}f Q dt'B{t')dt' 



(3) 



Typically the fluctuations of the magnetic field will occur in a time scale that is too short to be resolved by any 
measuring apparatus and thus one will have to resort to a time averaged description [13]. The fluctuations in the 
phase will induce a decay in certain coherences of the time-averaged density matrix. The evolution of the state vector 
will be as follows \tp{t)) = a |00) + /?e^W |01) + e i *(*) 7 1 10) + Se 12 ^ |11), but this is just for one realization and one 
needs consider the average density operator which would be given by: 

( \a\ 2 e -^Ma< 7 * e =2l W)a5* \ 

e^W(3a* \f3\ 2 /3 7 * e-W fiS* 

e #(*) 7 a* 7 /3* | 7 | 2 e-WfS* 



\ e 2 ^(*)<5a* eWSP* e^S-y* \S\ 2 J 

where the overbar denotes an ensemble average (and the ergodic hypothesis has been invoked to switch from a time 
to an ensemble average). It is reasonable to assume that <f> is normally distributed and thus e™ 41 ^*) = e l "^(*)-™ 2 <K*) 2 / 2 . 
The magnetic field B(t) can be modelled as white noise with zero mean which implies that <p(t) — and <fr(t) 2 = t/T 2 
where T 2 characterizes the variance of the random process B(t). Under this assumption it is clear that all coherences 
except the one proportional to /3* 7 approach zero as t ^S> T 2 , and the resultant state is an X state, with a = \a\ 2 , 
b = \/3\ 2 , c = | 7 | 2 , d = |(5| 2 , a = H 2 , z = and w = 0. 

The above example illustrates one of many possible ways in which an initially arbitrary pure state decoheres into 
X states. In the following section we intend to investigate general properties of these states. First in section (II) we 
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introduce and relate two convenient parameterizations for X states. In section (III) we calculate several important 
measures of quantum correlations and comment on the type of X states that maximize and minimize such correlations. 
We also include a measure of classical correlations for these states. Finally, in section (IV) we study the types of 
dynamics that generate and preserve the shape of X states. 



II. DEFINITIONS AND PARAMETERIZATIONS 



For a pair of two levels systems, A and B, one can define local bases with states that we label \0) A , \1) A for system 
A and |0) B , |1) B for system B. By taking the tensor product of the basis elements one can construct a basis for the 
qubit Hilbert space. In such a basis we write a general X state of the system as: 

px = a |00) (00| + b |01) (01| + c|10) (10| + d |11) (11| 

z |01) (10| + z* |10) (01| + w 1 00) (11| + w* |11) (00| . (4) 

In the ordered basis {|00) , |01) , 1 10) , |11)} the operator p x takes the matrix form shown in (1). Normalized density 
operators are positive semidefinite and have unit trace, which for the above parameterization implies the following 
constraints: 

a + b + c + d = 1, 
a, b,c,d> 0, 

\z\ < Vbc, and |tu| < Vad. (5) 

It is always possible to apply a pair of local unitary transformations to make all the coefficients in the definitions (4) 
and (1) non-negative. For instance, the phases of z and w can be absorbed into |0)^ and |0) B , i.e. one redefines 
e lar gO) |Q) A as |0a) and e 4ar s( M ') |0) B as |0s). Because the correlations of a system do not change when local unitaries 
are applied to the subsystems then such correlations will depend only on the absolute values of the coherences z and 
w but not on their phases. From now on we assume that such local unitary transformations have been applied and 
thus all the coefficients in (1) are real and non-negative. 

An equivalent and rather useful way of writing the operator (4) is by using the Fano parameterization[14]. To this 
end, one defines the usual one-qubit Pauli operators, 

a x = a x = |0) (1| + |1> <0| ; a 2 = a y = i (|1) (0| - |0) (1|) ; 
<7 3 =^ = |1)(1|-|0)(0|; <7 =I=|0)(0| + |1)(1|. 



With this definitions one can write: 

Px = 



^ 1 1 ® I + A 3 a 3 ® 1 + B 3 I® as + CiCTi ® 
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The two parameterizations can be related as follows: 

A 3 = (a + b)-(c+ d) B 3 = (a + c) - (6 + d) , , 

d = 2{z + w) C 2 - 2(z-w) C 3 = (a + d)-(b + c). W 

Notice also that from the Fano parameterization the reduced density matrices can be obtained straightforwardly: 

p A = tr B (px) = ^{l + A 3 a 3 }, p B =tr A (p x ) = ^{I + B 3 a 3 }. (7) 
Both density matrices are diagonal in the |0) , |1) basis. 

Finally, notice that X states contain as particular instances many important states like the four maximally entangled 
Bell states, 

|to> = -^(|00> + |ll», |&> = (I®*i)|tfo) (8) 
the maximally mixed state, mixtures of maximally entangled and maximally mixed states like Werner states 

Pw = (l-e)l+e\<f>i){<t>i\, (9) 



and all the states that for a given value of their mixedness (or purity) maximize their entanglement [8] . 
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III. QUANTUM AND CLASSICAL CORRELATIONS 



In this section we provide expressions for several quantum correlations for two-qubit X states. We also include 
an expression for the classical correlations of these states. Before calculating any relevant measure of correlations 
between the two parties in the bipartite system we write its state in two canonical forms. The first one is the eigenvalue 
decomposition p — J2i with: 

A 1/2 = u+± ^u 2 _+w 2 (10) 
= r+ ± ^Jr 2 _+z 2 

|Al/ 2 ) = ^r- ± y/u*_+W*} 100) + W |11)) 

|A 3/4 ) = ^({,- ±A A^}|01) + z|10>) 

where u± = (a ± c£)/2, r± = (6 ± c)/2 and iV]y 2 and Mjm arc normalization constants. Notice that (5) ensures 
that all the Xi are positive. From the above it is straightforward to obtain the von- Neumann entropy as S(px) = 
~Ei ^ lo g( A i)- 

Using the eigenvalues of px or directly calculating the trace of the square of the density matrix one obtains the purity 
of the state: 

trp^ = ^2 X 2 = a 2 + b 2 + c 2 + d 2 + 2w 2 + 2z 2 . 

i 

The purity will only be equal to one if either a = d = w = and be = z 2 or b = c = z = and ad = w 2 . 
From the definitions in this section, we calculate several relevant quantum correlations. 



A. Entanglement 

For decades one correlation, Quantum Entanglement, was considered the defining property that distinguishes quan- 
tum systems from classical ones. In simple terms, it was considered to be present in a system by the inability to factor 
its state into a product of the states of the subsystems that make it. To be more specific in quantifying entanglement, 
in the case of pure states of two qubits, this is achieved by calculating the von Neumann entropy of the reduced 
density matrix of the system. For mixed states, this method cannot be applied directly. First, the mixed state is 
written in a pure-state decomposition as follows: 

P = I>hfc>M- (ii) 

i 

Entanglement for this mixed state is then defined in terms of the entanglement of the pure states involved in the 
decomposition and minimized over all decompositions (since decomposition (11) is by no means unique): 

E{ P ) = minY J PiE{^i), (12) 

i 

Standing on this concept, Wootters defined the famous entanglement measure known as concurrence [15] . In the case 
of our X states, which in the most general cases are mixed, concurrence is given by: 

C(px) = 2 max jo, z — Vad, w — Vbc^ (13) 

Note that this quantity takes values between and 1, the former corresponding to no entanglement in the system, 
the latter to maximal entanglement, and anything in between correspond to cases of partially entangled states. 
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B. Partial Transpose and Negativity 



One of the most powerful techniques for entanglement detection is the use of positive but not completely positive 
(PNCP) maps. A positive linear map A between the space of operators acting in two Hilbcrt spaces Ha and Hb 
satisfies[16, 17]: 

If A{X) = Y then A(A f ) = F f 
If X > then A(X) =Y>0 

A positive map A is completely positive if for an arbitrary Hilbert space He the map Xq A is positive, otherwise is 
termed positive but not completely positive. Notice that a PNCP map will always map separable density operators to 
separable density operators. A failure to obtain a positive operator after applying a PNCP over a density operator will 
imply that the density operator that is given is entangled. One very useful example of a PNCP map is the transpose. 
In [18] the partial transposition of a density operator with respect to one of its subsystems (which is a PNCP) is used 
to give a very powerful condition for entanglement detection. In [19] it is shown that this condition is necessary and 
sufficient for entanglement in qubit-qubit and qubit-qutrit systems. Finally, in [20], this concept is used to quantify 
entanglement. A partial transpose is performed over the first subsystem (system A here), and the measure called 
Negativity, N(p), is the sum of the absolute values of the negative eigenvalues of the partially transposed matrix. For 
arbitrary dimensions if N{p) > there is entanglement in the system, but N{p) = gives us no definite answer as to 
whether the system in entangled or separable in the general case. Nevertheless, as already mentioned for qubit-qubit 
and qubit-qutrit systems, N{p) = necessarily implies that p is separable. The partial transposition amounts to 
swapping the second set of labels in each bra and ket in equation (4) which for X states amounts to exchange z and 
w. Thus the conditions for having a positive partially transposed X state are given [21]: 



z < Vad and w < Vbc (14) 

which is precisely what the concurrence tells us (13) since if any of these conditions is violated the partially transposed 
X state will not be positive or equivalently the state itself will be entangled. As for the precise eigenvalues of the 
partial transpose of an X states these are given by equation (10) but again doing the change w < — ► z. In this case 
only the cases with minus signs before the inequality might give negative eigenvalues and thus the negativity is given 
by: 



Af(p) = - min jo, u + - \J u 2 _ + z 2 , r + - yV 2 + u; 2 j 



where u± and r± are the same quantities that appear in (10). 

Another interesting property of a partially transposed qubit-qubit density matrix is that it will have at most one 
negative eigenvalue, and thus its determinant gives a necessary and sufficient condition for entanglement. In [22] this 
is shown with complete generality, here we show it for X states by elementary means. Notice that if for instance 
z > ^/ad then from (5) one has Vbc > z and \fad > w which automatically implies that y/bc > w. The argument is 
reversed if w > Vbc. 

C. Fully Entangled Fraction 

The fully entangled fraction (FEF) is defined as the maximum fidelity that a given quantum state has with a 
maximally entangled state [23]. The fidelity between a (generally) mixed state and a pure state is defined as [24]: 

jr(p,|v>>^|) = <V>|p|V> (15) 

For the fully entangled fraction one calculates the above expectation value with the four Bell states (8). For p x one 
has: 

£{p x ) = max(a + d + 2w - 1, b + c + 2z - 1). (16) 
For a state with a fixed value of the FEF £ it is easily found that the concurrence of such state is bounded by: 

£ < c < ^ (17) 
In particular the states that saturate the upper and lower bounds of the inequality are of the X type[25]. 
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D. The Schmidt Number of the State 

An important characteristic of a bipartite quantum state is its ability to encode information about a local quantum 
process that acts on only one of the qubits. In particular, one asks the question of whether by knowing an initial 
state pin and the state that is obtained by applying an unknown quantum process in one (say the second) qubit, 
Pout = (I <8> £) {An}, is it possible to know what the process C was? In [26] it is shown that the necessary and 
sufficient condition to provide complete information about the local process C is that the Schmidt number of the 
density matrix equals the square of the dimensionality of the system on which the process acts, in our case 2 2 = 4. 

The Schmidt number, which is a familiar concept for pure states, can be found for bi-partite density operators 
as follows. For any state p one can define the matrix T^ v = Tr{p(a fl ® cr„)}/2[27], where \x and v take the values 
(0, 1, 2,3), where cto is the identity. We may use the singular value decomposition (SVD) to rewrite in the form 

Si s ifj^9v \ where is a set of non-negative real numbers and ffft and g„ \ (i — 1, 2, 3, 4) are two sets of orthonormal 
vectors (i.e. X^/mVa^ = X^S^'s^ = Defining the operators F t = X^/m^ and G » = E^W, we 

find p = J2i s iFi ® which is the Schmidt decomposition of p. The operators Fi and Gi are orthonormal, i.e., 
tr (F^ Fj) = tr(GjGj) = Si j (they arc not, however positive, hence they do not represent states). The number of 
non-zero s t is the Schmidt Number of p. For X states, one finds the singular values are related to the coefficients 
appearing in the Fano decomposition as follows: 

Ci 

31 = ~2~' 

|C 2 | 
52 = "2"' 

y/\ + Al + Bl + Cl±D 

53/4 = 271 ' (8) 



where D = y] {1 + A% + B% + Cf ) 2 - 4 (C 3 - A 3 B 3 f. 



E. Quantum Discord and Classical Correlations 

As described in an earlier subsection, the quantification of entanglement is done using the von Neumann entropy. 
In short, the idea is to use the fact that randomness is introduced into the system when a quantumly correlated 
particle is ignored (by taking the partial trace over it). Another method to quantify the strength of the quantum 
correlations of the system is to use the difference between the effect of measurement on a classical system compared 
to a quantum system; i.e., to use the fact that measurements disturb quantum systems, but not classical ones. This 
is the idea that Ollivier and Zurek's Quantum Discord is based on [28]. Like Wootters concurrence, quantum discord 
is defined for two-qubit systems. Labeling the two subsystems by A and B, a measurement is performed on B. If 
a disturbance is detected, then that indicates the existence of quantum correlations, otherwise it implies that they 
are absent. The disturbance is quantified by using the mutual information function, which gives an indication of how 
much information is shared between A and B. The difference between the mutual information function before the 
measurement and after measurement defines the discord. However, the set of projectors that are applied on B have 
to be chosen so that they give the maximal value for the measurement-induced mutual information function. Notice 
that there is in general a more complex hierarchy of quantum correlations' quantifiers in which different types of 
measurement schemes are applied of which the Discord is a particular instance [29]. Although for concurrence, there 
is an analytic expression for X states, it is not possible to find an analytic expression for discord for the general X 
state [30]. In general the problem of the calculation of the discord can be cast into the solution of two transcendental 
equations as it is shown in [7]. However, Luo [31] was able to find one subclass of X states for which an analytical 
expression can be given. This class has maximally mixed marginals (MMM): 

PA = PB = \l (19) 
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MMM) = 



which implies from equation (7) that A3 = By = 0. For the MMM states the discord is given by [31]: 

\ {(1 - Id) - |C 2 | - |C 3 |) log(l - |d| - |C 2 | 

(1 - Id] + |C 2 | + |C 3 |)Iog(l - jCij + |C 2 | - 
(1 + |dj - |C 2 | + |C 3 |) log(l + |Cij - |C 2 | 
(1 + |d| + |C 2 | - |C 3 |)Iog(l + Idl + |C 2 | 
i-c, _ 1 + C. 



-|c 3 | 

|C 3 |) 
|C 3 |) 
|C 3 |)} 



log(l - C) - 



■ log(l + C) 



(20) 



with C = max{|Ci|, |C 2 |, |C 3 |}. Although, obtaining an analytical expression for the discord has been shown to be 
impossible for states more complex than the MMM states a great deal of advance can be achieved by characterizing 
the set of states that have zero discord. This was done in [32]. Following the convention used in [32], we assume that 
measurements are made on subsystem A instead of subsystem B. Then it can be shown that the set of states that 
have zero discord is given by [32]: 



Pcl =^2?k \i>k) (V'fel 



(21) 



where {|"0fc)} is an orthonormal basis set for subsystem A. Since the set ft is known, it is possible to define a geometric 
measure of discord by simply measuring the distance (square norm in the Hilbcrt-Schmidt space) of a given state to 
the closest state in the set Q. This was done in [32] and the result for X states is simply: 



Vf ] { P x) = - A min {Cl + Cl C\ + C\ + A 2 } 



1 



(22) 



= - mm 



{4 (w 2 + z 2 ) , (a - cf + (b - df + 2(w + z) 2 } 



The geometric discord when B is measured is simply obtained by replacing Ay by By in the first equality or swapping 
b and c in the second one. Finally, we also point out that an X state is non-discordant if and only if is fully diagonal, 
which is equivalent to saying that it only has two non-zero singular values (See (18)). 

Although for general X states, there does not exist an analytic formula for discord, it has been shown that there 
exists a set of projectors that will give accurate results [30]. These projectors are labelled by the term maximal- 
correlation-direction measurement (MCDM). Moreover, in [33], an expression for discord in the case of b = c is 
derived. Here, we derive a very similar result, but for the slightly more general case in which b and c are not 
necessarily equal. Using a related approach to these studies, we find an expression that is a very good approximation 
to discord: 



Q(px) « S(p B ) ~ S( PX ) + min{JVi, N 2 } , 



(23) 



where S(px) is the von Neumann entropy of the general X state density matrix, S(pb) is the von Neumann entropy 
of the reduced density matrix of the second qubit (labelled B, on which the measurement is made) and, 



Ni=H 



No 



2 

alog 2 



l + y(a-d + b-c) 2 + 4(z + wy 



a 



a + c 



- b\og 2 



b + d 



- clog 2 



a + c 



dlog 2 



b + d 



(24) 



where H{y) = —y log 2 (y) — (1 — y) log 2 (l — y) is the binary entropy function. In [34], a parameterization, in 
terms of variables 9 and <f>, that is used to find discord is defined. In the language of this reference, Ni is found for 
9 = tt/4, </> = 0, and N2 is found for 9 = 7r/2, <j> = 0. To elaborate on the accuracy of this approximate expression for 
Q, when we analyze the 1 x 10 5 randomly generated X states (see [34] for method description), we find the following: 
% of the points have error > 10~ 3 , 0.001 % have error > 10~ 4 , 31.44 % have error > 10~ 5 , 86.10 % have error 
> 10~ 6 , and 90.78 % have error > 10~ 7 . This shows that if we are only interested in the accuracy up to the forth 
decimal place, then this is a very good approximation for discord of an X state. 

Notice that equation (23) also gives an approximate expression for the amount of classical correlations that a 
given state has. Recall that the quantum discord as defined by Zurek and Ollivier is Q(p) = I(p) — C(p), where 
I(p) = S(pa) + S(pb) — S(p) is the mutual information function, and C(p) is the measurement-induced mutual 
information function maximized over all measurements on subsystem B. One can also interpret these quantities as 
follows: I(p) represents the total (classical and quantum) correlations present in the system, Q(p) represents the 
quantum correlations, and C(p) represents the classical correlations (See [31, 35], for example). In the light of this 
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interpretation, we can use C(p) that we calculated to obtain the Q in (23) to represent the classical correlations in 
the X states. It is approximately given by: 

C(p x )^S( PA )-mm{N 1 ,N 2 }, (25) 

where S(pa) is the von Neumann entropy of the reduced density matrix of subsystem A, and Aq and N 2 are defined 
in (24). 

F. Measurement-Induced Disturbance 

Quantum Discord involves finding a local projector (on subsystem B) that minimizes its value. Moreover, Discord 
is not symmetric, meaning that the discord of A with respect to B is not the same as that of B with respect to A. 
That is why one has to think which subsystem is more appropriate to perform the local measurement on in order to 
define discord. These deficiencies inspired Luo to introduce the Measurement-Induced Disturbance or simply MID 
[36]. 

Like Discord, MID exploits the fact that measurements of quantum systems disturb them in order to capture 
quantum correlations. However, there are two differences. Now, the local measurements are performed on both 
subsystems A and B. In addition, MID does not require searching for the optimal set of local projectors. Instead, the 
chosen projectors are constructed using the eigenvectors of the reduced density matrices of A and B. MID is defined 
as follows [36, 37]: 

MID = I(p) - I(P(p)), (26) 

where 



p (p) = E E ® n f ) p ( n ^ ® n f ) - ( 2? ) 

»=i j=i 

and Uf and H? are projectors constructed using the eigenvectors of the reduced density matrices of systems A and 
B, respectively. Note that for two-qubit systems, and hence the X states we are focusing on in this paper, m, n = 2. 
In fact, for X states, P(p) is given by: 



Px{ P ) 

and the MID is given by: 



/ a \ 

6 

c 

Vo d) 



(28) 



MID X = -S(px) + S(P x (p)), (29) 

where S(px) is the von Neumann entropy of the A-matrix before measurement, and Px{p) is the von Neumann 
entropy of the density matrix after the measurement as given in (28). 

Note, however, that although MID is easier to calculate than discord, for several cases in which there is no quantum 
advantage, it predicts maximal quantum correlations [38, 39]. The approach taken in [39], in which Ameliorated MID 
is introduced, addresses this issue. For the X states, it is also shown that MID is ambiguous for the states classified 
as MMM. However, a simple calculation reveals that if Ilf and are constructed using the eigenvectors of a z , then 
MID — Q, where Q is given by eq.(23). It is also shown in [39] that for these states, MID and discord are the same 
for the Werner states and for the pure states. This implies that the result in (29) can be used in these two cases 
without ovcrcstimation of the strength of quantum correlations. Another related observation[40] in a class of states, 
also named Werner states, is the equality between discord and MID. Note, however, that the latter are not the same 
as the Werner states we discuss in this work, which are states that are a combination of the identity and a Bell state 
(See (9)). In [40] , the Werner states are defined to be those that satisfy: p = U®U pW ®W for any unitary operator 
U. Nevertheless, states in the two different classes have equal discord and MID. 
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IV. DYNAMICS 

In this section we study the types of dynamics and quantum channels that preserve the shape of an X state. One 
of the reasons why X states became so popular in the study of the dynamics of quantum correlations is because of 
the relatively simple form that the concurrence takes for such states. Moreover since in such studies [41-43] they stay 
in the X form for all times, the simple equation (13) is valid throughout their evolution. To study what types of 
dynamics preserve the shape of X states we employ the very elegant algebraic characterization of X states presented 
in [44] . The key ingredient of the characterization is to notice that the set of operators 

S ={I ® I, <7 3 ® I, I ® 03,(71 ® <7l,<72 ® (72, O3 ® &3, &2 ® & 3 , (7 3 ® (7 2 } (30) 

is closed under multiplication, i.e., the product of two of them will be proportional to another element of the set S. 
The other 8 operators that complete the 16 element basis for the 2 qubit Hilbert-Schmidt space are simply: 

S' ={I ® <7l,I ® (7 2 , CTl ® I, £72 ® I, <7l ® (7 3 , (7 2 ® (7 3 , CT 3 ® (7 2 , CT 3 ® CTl} (31) 

The product of operators belonging to S and S' satisfy the following properties: Let A, B e S and C,D £ S' then 
AC G S', AB e S and CD e S. With this elementary observations in mind we characterize the types of unitary 
evolution and non-unitary evolution that map X states to X states. 

For the Hamiltonian dynamics given by the von-Neumann equation, J^p = i[p } H) it is easily shown that p will remain 
X shaped if H is also X shaped. Equivalently p will remain X shaped if and only if H is spanned by S. For 
non-unitary evolution we first study the case of continuous time evolution. In this case the dynamics is given by a 
Markovian master equation in the Lindblad form [45]: 

d n 2 -i 

—p = i[p,H}+ h n,m {2L n pL\ n - pLl n L n - L\ n L n p) , 

n,m— 1 

where N is the dimensionality of the Hilbert space in which p acts, N — 4 in our case. H dictates the Hamiltonian 
dynamics that was studied in the previous paragraph. L n are a set of orthonormal operators in the Hilbert-Schmidt 
space to which p belongs. Notice that there are ./V 2 of such operators, but because one of them can always be chosen 
to be the identity (which will not cause non- unitary dynamics) the sum can be restricted to iV 2 — 1. Finally, h n , m 
are the complex entries of a positive semi-definite matrix. In the non-unitary part of the master equation the density 
operator appears multiplied by two different operators, L n and L m . These operators will preserve the X shape if and 
only if L n ,L m <G span{5} or L n ,L m e span{5'}. If cither the operators L m ,L n contain elements of both S and S' 
or they are spanned by different sets (S and S') then there will be products involving two operators from S and one 
from S' which, in general, will be outside the space of X states. Finally, as for quantum channels these are defined in 
the operator sum representation as [3]: 

p^p> = C{p} = Y,X iP X\ (32) 

i 

with XiXj = I. In this case, the same argument applied to the master equation can be used, that is, £ will map 
X states to X states if and only if the Xi do not mix operators from S and S' . 

V. CONCLUSIONS 

In this paper, we aimed to present some of the most interesting and useful quantum properties in the literature, 
calculated for X states. This work served as a reminder of already known results for these states, such as those 
pertaining to concurrence. We also presented some new results. First, by uncovering the singular values for these 
states, we have provided a tool to determine whether they can be used in ancilla-assisted state tomography. We also 
calculated measures of entanglement, other than concurrence, which can be of interest. This includes the Negativity 
and the Fully Entangled Fraction. Moreover, we derived results about the quantum discord of these states, one with 
regards to geometric discord and the other an approximate analytic expression for discord. The latter is shown to 
be very accurate when the last digit to be considered is the fourth one after the decimal place. Since it was shown 
that an analytic expression for discord of states that fall under the class of X states does not exist, this approximate 
(analytic) result for discord provides an easy and reliable tool in situations where optimization is neither practical nor 
necessary. Moreover, the measurement-induced disturbance, is calculated for X states, with the added caution that it 
can only be used for certain states such as the Werner and pure states, and that for the MMM states, the projectors 
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have to be chosen carefully. The X state MID turns out to be a simple expression, dependent on the entropy of the 
total system before and that after the measurement is performed. Finally, we discussed the dynamics that preserves 
the form of the X states by providing an exhaustive classification of the types of quantum master equations and 
quantum channels that preserve the shape of an initial X state. 
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